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Anisotropic antiferromagnets in an external magnetic field show a rich variety of different ground 
states meeting in transition lines and multicritical points. We study the dependence of the ground 
states of these systems in the three dimensional space on physical parameters as exchange, single 
ion and cubic anisotropy. 

One identifies four different ground states: the paramagnetic (PM), the antiferromagnetic (AF), 
the spin flop (SF) and the biconical (BC) ground state. In the case of absence of a cubic anisotropy 
the transition lines separating the different ground states can be calculated analytically, otherwise 
they have to be calculated numerically. We also considered the behavior of the staggered magneti- 
zation which characterizes the different ground states. From its behavior the order of the transition 
from one state to the other is determined. 

But also the order of the transition changes along the transition lines when including the cubic 
anisotropy, especially at the reeentrant region where a transition from SF to BC and back to SF 
by increasing the external field H occurs. Multicritical points are founded which are assumed to be 
tricritical or critical endpoints. 

The results obtained may be relevant for other systems since the antiferromagnetic model can be 
mapped to a lattice gas model where the biconical ground state is interpreted as supersolid phase. 
Recent renormalization group calculations show that such a phase would indicate the existence of a 
tetracritical point. 



INTRODUCTION 

Anisotropic systems in an external field not coupled to the order parameters (OP) characterizing the condensed 
phases show interesting phase diagrams depending on the values of the anisotropy and the external field. One example 
is the anisotropic antiferromagnet in an external magnetic field [1]. Other ones are crystal systems undergoing a 
displacive phase transition under external pressure [2]. The possible variety of phases at finite temperature is also 
determined by the topology of the ground state phase diagram since the phases corresponding to the different ground 
states extend to finite temperatures. However at finite temperature new phases not present at T = may arise [3]. Due 
to the variety of possible ground states or phases multicritical points of different kind are possible at finite T, Using 
the magnetic language this systems show antiferromagnetic (AF), spinflop (SF), biconical (BC) and paramagnetic 
(PM) ground states and phases. The transition lines between these phases may meet in one point: (i) the bicritical 
point, where three phases coexist (AF, SP, and PM) and (ii) the tetracritical point, where four phases (AF, BC, SP, 
and PM) coexist. Going beyond mean field theory [3] the nature of the multicritical points has been studied within 
renormalization group theory [4] (in a one loop order calculation) . These results have been questioned by higher loop 
calculations [5, 6]. Recently both using simulation methods [7, 8] and renormalization group theory [9, 10] the nature 
of the multicritical point have been reconsidered. 

It is well known that there is a correspondence between such magnetic model and the quantum lattice gas [11] which 
allows to transpose the results obtained for the magnetic system to systems like He 4 with superfluid and supersolid 
phases [12]. The question of the existence of such phases is of general interest. 

Symmetry consideration are essential for the existence and the nature of ground states and phases. This is also the 
case for the systems treated here. It has been found that interaction terms of cubic symmetry has strong influence 
on the phase diagram [13]. Here we present results at T = which clarifies the order of the transition lines between 
the different ground states. Moreover we extend the models studied so far by combining several types of anisotropy. 
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Table I: Correspondence between the phases of the classical magnetic system and those of the quantum lattice gas (QLG). Also 
shown are the order parameters, the alternating magnetization in the fled direction iVy and perpendicular to it N±, and the 
components of the magnetization My and M± 
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SINGLE ION ANISOTROPY 



Let us first consider the well studied case of an antiferromagnet with the Hamiltonian H 

N A ,N B I N A N B \ I N A N B \ 

n = j ( A + S iy S jy ) + S iz S jz ) - h £ s« + E s i* + D E So + E s%) (i) 

It describes Na and Nb interacting classical spins «S, on two sublattices A and B, where Na = Nb- The anisotropy 
of the positive exchange interaction J is characterised by the parameter A with A = for a magnet where only the 
^-components of the three dimensional spins interact and A = 1 for an isotropic magnet. The spins are subjected to 
an external magnetic field H and a single- ion anisotropy D. 

The case where beside the exchange anisotropy only a single-ion anisotropy is present can be solved analytically [14] 
at T — 0. ABC ground state exists for a certain region of positive D and external magnetic field H for A = 0. This 
region becomes smaller for increasing A and diminishes for an isotropic antiferromagnet with A = 1. The transition 
lines between the AF to the PM and between the AF and the SF ground state are of first order. The other transitions 
are of second order. Figure 1 shows the whole ground state diagram with the parameters chosen so that a BC ground 
state appears. Here we divide the parameters H,D,F by the number of next neighbors (2 = 6 for simple cubic 
system), so in order to compare with the results in [13, 15] one has to multiply all the three parameters in this paper 
with 6. Also the OP for the complicated BC ground state can be calculated exactly which reads: 
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An extension to finite temperature within mean field theory has been discussed in [3]. As already explained in [11] 
and used in [14] a quanta! lattice gas can be mapped on an anisotropic antiferromagnet with a Hamiltonian of form 
(1). The correspondence of the different ground states is given in Tab. I. 

In the H — T ground state diagram the ground states at T = extend as phases to higher temperature. Since for 
high enough temperature and large enough external magnetic field the PM phase should be the stable one, all other 
phases are confined to finite regions in the phase diagram. Thus it happens that the phase transition lines cross each 
other in multicritical points. This might be a bicritical point when the AF, the SF and the PM phase meet, or a 
tetracritical point where the AF, the SF, the BC and the PM phase meet. It should be noted that this is even possible 
if the BC phase does not extend to T = [3]. 



3 



D=0, A=0.8 



D=0, A=0.8 



I 1.0 




- 1 a 



1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 
F 



Figure 2: (online color) Ground states minimizing the energy, Eq. (3), in the magnetic field H and cubic anisotropy F plane 
at an exchange anisotropy A = 0.8 and with zero single ion anisotropy D. The ground states are denoted as follows: AF 
antiferromagnetic, PM paramagnetic, SF spin flop ; BC1,2 biconical oriented along the diagonal or the axis in the x — y-plane. 
Shown are special points where either more than two ground states meet or where the order of the transition line changes from 
first (dashed line) to second order (solid line). These are the triple points TR1,2, the tricritical points Tl,2 and the critical end 
points CE1,2. 



CUBIC ANISOTROPY 



Uniaxial antiferromagnets may also have cubic anisotropy. Here we consider an anisotropic antiferromagnet with a 
cubic anisotropy F instead of the single-ion anisotropy D. For such a system the Hamiltonian reads 

Na-Nb I N A N B 

H = J ^ ' (A (Si X Sj X + SiySjy) + SizSjz) H I S ' S% x + S ' S j z 



Na N b 

- F ^2 ( S ix + S iy + St*) + X! ( S % + S Jv + Sjz) 



(3) 



The effect of such a type of the cubic anisotropy has been studied recently by [13, 15]. It has been shown that two 
types of BC ground states are possible (BC1 for F < and BC2 for F > 0), depending whether the perpendicular 
component of the alternating magnetization lies in the diagonals or along the axes in the xy-plane. Moreover it has 
been found that most of the transition lines between the ground states become first order. There was also found a 
reentrant region for the SF ground state becoming unstable for increasing the external magnetic field with respect to 
the BC state, which for even larger field again becomes unstable with respect to the SF phase (see Fig. 2). In this 
reentrant region the transition line has been found to be of second order. 

We verify the occurrence of the reentrant region [13] by calculating the order parameter of the ground states which 
take place in the vicinity of this region (see Fig 3). The 2D figure shows the the SF-BC1 occurs by increasing the 
applying field H. Continue to increasing the field further the SF reappears again. 

In order to find out the transition points where a change in the order of the transition takes place we have to 
calculate the nonzero order parameter component along the border lines for the case treated in [13]. Along the 
transition line between the BC1 and the PM ground state the parallel staggered magnetization N\\ goes to zero at the 
point where the SF - PM and the BC1- SF transition lines meet (see Fig. 2). The behavior of the OP is compatible 
with a power law with a correction term 
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where F* is the value of the cubic anisotropy at the special point, and F(H) the value of the cubic anisotropy along 
the transition line. A, B and F* are fit parameters. This point is a critical end point (CE1). From this critical 
endpoint the second order reentrant line starts. This line ends in a point where the AF, the BC1 and the SF ground 
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Figure 3: (online color)A 3D and a 2D diagram showing the occurrence of the reentrant by increasing the magnetic field H 
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Figure 4: (online color)Left: Detail from the right part of Fig. 1 on a larger scale around the reentrant region. Solid lines are 
of second order, dashed lines of first order. The second order reentrant line between the biconical and spin flop ground state 
starting at the tricritical point TR1 ends in the critical end point CE1 on the transition line to the paramagnetic ground state. 
On the transition line between the paramagnetic and spin flop ground state one finds the tricritical point Tl. Right: Comparison 
of the transition line between the SF and PM ground state calculated on the assumption of a second order transition (straight 
line) and the transition line calculated numerically by minimizing the energy. Indeed when the transition curve deviates from 
the second order line the transition becomes first order. 



state meet. This point turns out to be a triple point (TR1) since the second order character of the reentrant line is 
changed to first order (at T2) before it reaches the triple point. The value of Nn becomes nonzero on the reentrant 
line at T2 and this point might be a tricritical point (see Fig. 5). 

We also find a change of order along the transitions line between the SF and the PM ground state. At Tl the 
transition becomes second order for larger values of F(H). This is just at the point where the transition line begins 
to deviate from the analytically calculated second order transition line (see Fig. 5). 

Conclusively we summarize all the critical points in a table (see Tab II) 



Table II: (online color)Summary of all the special points discovered in the ground state diagram for D = and A = 0.8. There 
exist tricritical points (T), tripel points (TR) and the critical end points (CE). 
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(0,0.6) 
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Figure 5: (online color)Left: Change of the parallel order parameter |7V|| | going along the first order transition lines between 
the BC1 and PM ground states. The dots are calculated values and the solid line a fit according to Eq (4). The zero value of 
the parallel OP is reached at the critical end point CE1. Right: Similar the change of the perpendicular order parameter |JVj_| 
going along the first order transition lines between the the SF and PM ground state. The zero value of the perpendicular OP 
is reached at the tricritical point Tl. 



SINGLE ION AND CUBIC ANISOTROPY 



Let us now consider the general case when the exchange, single-ion and cubic anisotropy are present 
H 



N A ,N B I N A N b 
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(5) 



Then the planar ground state regions extend to three dimensional regions and whole lines of multicritical points 
appear. In Fig. 7 we give as example the case where the single ion anisotropy has been chosen nonzero (A = 0.8) 
and a sketch of the three dimensional ground state diagram. An analytic result for the transition line between the SF 
ground state and the PM ground state can be calculated assuming a second order transition, which of course is just 
correct in certain section of the transition line (see Fig. 4). The relation reads 



H(A, D,F) = 2D + 4F + A + 1 



(6) 



Also in zero magnetic field (H = 0) the stable ground states as function of the anisotropies can be calculated 
analytically. The transition between the AF and the BC2 ground state starts at a positive value of the cubic anisotropy 
F 



F(D, A) = (1 - 2D- A)/4. 



(7) 



The transition occurs just only in the region D < i-^A and is continuous (see Fig 6), where the OP parameter of the 
biconical ground state (BC2) can be given. 

8a b — arccos ( ± 



1 + 2F - 2D 



6F 



(8) 



This 2nd order transition continues to be valid also when applying when magnetic field. Increasing the single ion 



isotropy ( D > i-^A ) the Spin-Flop is preferred to the BC2 and the transition takes place at 



F(D,A) 



-1 + 2D 



(9) 



If the cubic anisotropy is zero no transition to the BC1 ground state takes place. 
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Figure 6: (online color)Left: The analytically calculated ground state diagram for H = 0. Right: The continuous transition 
from AF to BC2. The two polar angles, from which the OP can be calculated, of the two different sublattices are shown 
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Figure 7: (online color)Left: Ground states minimizing the energy, Eq. (7), in the magnetic field H and cubic anisotropy F 
plane at an exchange anisotropy A = 0.8 and with single ion anisotropy D = 0.05. Right: A 3D-phase diagram H (F, D) with 
fix single-ion-anisotropy A = 0.8 

For negative cubic anisotropy in zero magnetic field the transition from the AF to the SF state takes place at 

D(F,A) = i(l-A) (10) 

and is independent of F. 

CONCLUSION AND OUTLOOK 

A system of an anisotropic antiferromagnet shows lots of interesting ground states and critical points. It can also 
mapped into a quantum lattice system where the biconical ground state is interpreted as the supersolid state. We 
have shown that for the simple case where just the exchange (A) and the single-ion anisotropy (D) are considered one 
can calculate analytically all the ground states. In the case with cubic anisotropy the resulting Hamiltonian is more 
difficult to be solved, so that we have to calculate the ground states numerically. In this system two kind of biconical 
ground states appear, the BC1 for F < and BC2 for F > 0. 

Most of the transition lines are of first order, but in some cases the order of the transition also changes along 
the transition line. So lots of interesting critical points appear ( critical endpoint, tricritical point). E.g along the 
transition line between BC1 and SF for F < a reentrance region can be shown where the order of the transition 
changes from second to first order by increasing the absolute value of the cubic anisotropy. 
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However in the absence of the field H we can derive the transition surfaces analytically, because the term with the 
applied field (H) is the part of the Hamiltonion which breaks the symmetry and therefore causes this Hamiltonian 
unsolvable. The BC2 ground state does exist also at H = for large value of F. 

The topology of the ground state diagram at T = determines the possible multicritical points at finite temper- 
atures. Without cubic anisotropy bicritical or tetracritical points are possible where two lines of second order phase 
transitions with order parameters of dimension n = 1 and n = 2 meet. Cubic anisotropy may change this picture and 
multicritical point where two lines with n = 1 and n = 3 or n = 3 and n = 2 might meet (see Figs. 2,7 for positive 
values of F). Cubic anisotropy changes most of the transitions to first order, however on the transition line separating 
the BC1 from the SF ground state part of this line starting at a critical end point is of second order type. Moreover 
part of the transition line separating the spinflop state from the paramagnetic state is of second order. 

We thank W. Selke for valuable discussions. This work was supported by the Fonds zur Forderung der wis- 
senschaftlichen Forschung under Project No. P19583-N20. 
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